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Technical  Report  Summary 

This  is  the  fourth  semiannual  report  dealing  with  an 
investigation  of  multiple  seismic  events  and  first  zone 
discriminants.  Reported  here  are  the  results  obtained 
from  studying  the  ratio  of  amplitude  spectra  and  cross- 
correlation studies  in  the  time  domain. 

Based  on  first  zone  measurements  and  the  application 
of  the  appropriate  equations,  calculated  source  radii  were 
found  to  be  clearly  different  for  the  different  type  sources. 
On  the  basis  of  this  radius  study,  the  source  dimension  of 
an  earthquake  is  significantly  larger  than  the  source 
dimension  of  a nuclear  explosion,  for  the  same  order  of 
magnitude  m^.  It  should  be  noted,  however,  that  the  dif- 
ference in  source  radii  results  from  application  of  the. 
correct  equations  and  does  not  result  from  a difference 
in  peak  frequency.  Hence,  equivalent  source  radii  does 
not  appear  an  unambiguous  method  of  differentiation.  In 
addition,  this  kind  of  test  suffers  from  the  failure  of 
not  being  able  to  be  performed  by  automated  devices. 

The  slope  of  the  displacement  spectra  in  the  range  of 
frequencies  from  .20  Hz  to  the  corresponding  peak  frequency 
are  clearly  different  for  natural  events  and  nuclear 
explosions.  The  slope  for  this  range  of  frequencies  is 
near  zero  for  the  Massachusetts  Mountain  earthquake  and 
clearly  positive  for  all  unknown  nuclear  events  and  the 
double  detonation  BLENTON /THISTLE . The  fact  that  the 


slope  of  the  displacement  spectra  is  positive  is  very 
significant  when  compared  with  the  results  obtained  by 
Wyss  and  others  (1971)  for  the  single  detonation  MILROW, 
suggesting  that  the  multiplicity  of  BLENTON/THISTLE  appears 
to  have  no  influence  on  the  general  behavior  of  the  con- 
sidered part  of  the  displacement  spectra.  The  advantage 
of  this  potential  discriminating  method,  under  the  concept 
of  an  unmanned  observatory,  is  that  of  being  amenable  to 
analogue  methods. 

While  the  slope  of  the  displacement  spectra  for  the 
long  periods  appears  to  be  a good  method  of  identification, 
it  suffers  in  that  broad  band  seismic  data  is  required. 

This  limits  the  utility  of  uncorrected  data  from  the  rela- 
tively narrow  band-pass  short  period  Benioff  system. 

Attempts  to  correct  the  Benioff  data  used  in  this  study 
were  unsuccessful  for  frequencies  less  than  0.2  Hz. 

The  spectra  of  the  seismic  background  noise  in  the 
area  of  study  appear  to  have  a strong  peak  around  .16  Hz, 
as  has  been  established  by  Brune  and  Oliver  (1959)  and 
Brune  (1971). 

The  modulations  preser,t  in  the  amplitude  spectra  appear 
to  be  very  complex.  The  use  of  the  theoretic?!  modulation 
functions  were  successful  only  for  events  BLENTON/THISTLE 
for  KNB  and  LAC  stations,  A for  MNV  station,  B for  ELK, 

LAC,  and  MNV  stations,  C for  ELK  and  KNB  stations,  D for 
LAC  stations,  and  E for  LAC  station.  For  the  remaining 
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stations,  no  reasonable  results  were  obtained.  The  modula- 
tion functions  appear  to  be  highly  sensitive  to  the 
variation  of  the  delay  times.  The  use  of  more  involved 
theoretical  modulation  functions  will  increase  this 
sensi'ivity  and  therefore  be  of  doubtful  value  as  a sole 
approach  to  determine  the  multiplicity  of  a source,  but 
nevertheless,  can  complement  other  methods.  The  attempts 
made  to  determine  the  theoretical  envelopes  of  the  amplitude 
spectra  ratios  were  unsuccessful  due,  in  part,  to  the 
number  of  parameters  involved  and  the  strong  dependency  on 
frequency  of  these  theoretical  modulation  functions.  In 
addition,  for  the  purpose  of  this  study,  events  A through 
E were  assumed  to  have  an  unknown  character  (single  or 
multiple  c/ent).  If  they  were  multiple  events  it  was  not 
known  what  the  source  geometry  or  delay  times  were. 

Higher  values  of  the  correlation  at  zero  lag  were 
obtained  when  nuclear  events  were  correlated  with  nuclear 
events  than  when  a natural  event  was  correlated  with 
nuclear  explosions.  Although  the  value  of  .08  for  the 
correlation  of  the  single  explosion  DIDO  QUEEN  with  the 
double  event  BLENTON /THISTLE  is  anomalously  low,  conclu- 
sions should  not  be  derived  from  the  results  of  only  one 
station.  Although  an  exhaustive  study  was  not  attempted, 
it  appears  that  the  correlation  at  zero  lag  is  useful  as 
the  actual  correlation  waveform. 
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The  consistent  high  values  obtained  for  the  correlation 
of  Event  B with  the  double  detonation  BLENTON/THISTLE 
suggest,  on  the  basis  of  the  theoretical  frame  considered 
m this  investigation,  the  possibility  of  a multiple  source 
for  Event  B.  On  the  same  theoretical  basis,  the  consis- 
tently low  values  obtained  for  the  correlation  of  Events 
A and  C with  the  known  double  explosion  BLENTON/THISTLE 
indicate  the  possibility  of  a single  detonation  for  both 
events.  The  significant  differences  in  the  values  of  the 
correlation  at  zero  lag  of  Events  D and  E with  BLENTON/ 
THISTLE  ban  any  suggestion  about  the  nature  of  the  source. 

Although  the  results  obtained  with  the  correlation 
analysis  suggest  the  method  as  a potential  tool  for  the 
discrimination  of  nuclear  events  on  the  basis  of  the 
multiplicity  of  their  source,  it  suffers  from  the  diffi- 
culty of  not  being  adaptable  to  analogue  methods  due  to 
the  large  number  of  reference  signals  required. 
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Dear  Sir: 

This  report  is  a summary  of  research  dealing  with  multiple 
seismic  events  and  first  zone  discriminants.  The  research 
is  divided  into  the  following  categories  and  will  be  dis- 
cussed individually. 
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Introduction 

An  underground  nuclear  explosion  will  generate  seismic 
waves  that  can  be  detected  and  recorded  providing  the  yield 
is  large  enough  and  the  distance  between  source  and  seismo- 
graph is  not  too  large.  During  the  past  approximately  18 
years  considerable  effort  has  been  directed  toward  the 
development  of  seismological  methods  to  detect  underground 
nuclear  events  and  to  discriminate  between  earthquakes  and 
nuclear  explosions. 

For  seismic  events  at  teleseismic  distances,  Willis 
and  others  (1963)  approached  the  problem  in  terms  of 
amplitude  characteristics.  The  same  authors  (1963)  carried 
out  a study  of  the  spectral  differences  of  earthquakes  and 
nuclear  explosions.  Brune  and  others  (1963),  Press  and 
others  (1963) , Liebermann  and  others  (1966) , Evernden 
(1969)  and  Liebermann  and  Pomeroy  (1969)  have  investigated 
various  criteria  which  are  mainly  related  to  the  detection 
of  depth  and  excitation  of  surface  waves  recorded  on  long- 
period  seismographs.  On  the  basis  of  these  and  other 
studies,  it  is  generally  accepted  that  given  an  earthquake 
and  a nuclear  explosion  releasing  equivalent  amounts  of 
energy,  the  earthquake  will  generate  more  long-period 
surface  waves  than  the  nuclear  detonation.  For  teleseismic 
events  of  magnitude  above  4.75  and  possibly  £4  (Evernden 
and  others,  1971),  the  most  reliable  method  appears  to  be 
the  difference  in  the  average  ratio  of  surface-wave 
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magnitude  (Ms)  to  body-wave  magnitude  (m^) , which  in  general 
is  smaller  for  nuclear  explosions  than  for  earthquakes.  When 
this  identification  method  fails,  there  are  other  criteria 
that  could  be  applied,  which  are  based  on  an  accurate 
depth  determination,  the  different  complexity  of  the 
recorded  seismic  waves  and  the  ratio  of  P-waves  energy  to 
S-waves  energy.  Therefore,  at  the  present  time,  the 
identification  at  teleseismic  distances  of  single  seismic 
events  with  magnitudes  mb  4.75  appears  to  be  feasible 
with  a high  degree  of  reliability. 

When  the  energy  released  by  a seismic  event  is  less 
than  that  resulting  from  a magnitude  4.75  event,  the 
reliability  of  the  mentioned  methods  decreases.  In  addi- 
tion, the  situation  becomes  more  complex  when  multiple 
sources  are  considered.  Furthermore,  it  has  been  recog- 
nized that  the  degree  of  differentiation  that  can  be 
obtained  by  seismological  methods  is  strongly  dependent 
upon  spatial  distribution  and  density  of  detectors  and 
upon  the  distance  of  the  receiving  stations  from  the 
unknown  event.  The  net  result  is  that  at  the  present  time, 
identification  cannot  be  achieved  with  the  required  degree 
of  confidence.  In  full  awareness  of  this,  nations 
interested  in  a comprehensive  treaty  continue  research 
on  the  subject. 

The  majority  of  the  effort  to  date  has  involved  the 
identification  of  seismic  events  at  teleseismic  distances. 
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For  smaller  yield  events,  especially  small  multiple  source 
events,  it  is  necessary  to  resort  to  seismic  data  recorded 
at  first  zone  distances.  First  zone  distances  as  used  in 
this  report  are  distances  less  than  1000  km  between  the 
source  and  the  receiver . 

Notable  differences  in  the  amplitude  spectra  of  earth- 
quakes and  nuclear  detonations  in  the  frequency  range  from 
.02  Hz  to  0.5  Hz  have  been  reported  to  exist  by  Pasechnik 
in  the  International  Institute  for  Peace  and  Conflict 
Research  Report  of  August,  1968.  Wyss  and  others  (1971) 
conducted  an  investigation  of  the  P-waves  spectral  dif- 
ferences in  the  frequency  range  from  .03  Hz  to  2.0  Hz  for 
natural  events  and  underground  nuclear  explosions,  recorded 
at  teleseismic  distances  with  magnitudes  (m^)  from  5.6  to 
6.5  and  from  5.9  to  6.6,  respectively.  The  results  of  this 
investigation,  as  reported  by  the  authors,  indicated  a marked 
difference  in  the  behavior  of  the  displacement  spectra  of 
natural  and  explosive  sources  and  this  suggests  a potential 
discriminant  for  both  types  of  events.  It  does  not  appear, 
however,  that  these  spectral  differences  have  been  investi- 
gated at  first  zone  distances. 

Based  on  the  results  obtained  by  Wyss  and  others 
(1971) , a similar  study  will  be  conducted  using  eight 
nuclear  explosions,  one  of  them  known  to  be  multiple,  and 
one  earthquake  from  the  Nevada  Test  Site.  The  recording 
stations  utilized  will  all  lie  within  the  first  zone. 
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Since,  in  the  amplitude  spectra  of  a seismic  signal 
the  phase  formation  is  lost  and,  furthermore,  the  spectra 
is  dependent  upon  the  propagation  path  of  the  seismic 
energy,  cross-correlation  studies  in  the  time  domain  will 
also  be  conducted  in  conjunction  with  the  ratio  of  the 
amplitude  spectra. 

In  summary,  the  purpose  of  the  research  effort  reported 
here  is  to  investigate  the  potential  of  spectral  analysis 
and  correlation  techniques  as  first  zone  discriminants 
between  nuclear  explosions  and  earthquakes.  This  report 
is  based  on  the  thesis  prepared  by  Tarazona  (1975)  which 
was  sponsored  by  this  project. 
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Theory 

Amplitude  Spectra 

The  mechanisms  of  earthquakes  and  explosions  have 
been  studied  and  modeled  theoretically  by  several  investi- 
gators (Sharpe,  1942;  Kasahara,  1957;  Archambeau,  1964; 
Savage,  1966;  Berchkemer  and  Jacob,  1968;  Aki,  1968; 
Haskell,  1964,  1969;  Brune,  1970,  1971).  In  all  models 
with  a faulting  mechanism,  there  exists  a frequency,  fQ, 
known  as  the  corner  frequency,  which  is  related  to  the 
source  dimension  through: 


where  c is  a constant  whose  value  is  near  1 and  depends 
upon  the  specific  type  of  model  used;  v is  the  velocity  of 
the  body  waves  (P-waves  or  S-waves  velocity,  depending  on 
the  phase  used  in  the  far-f ield  spectrum  calculation) ; and 
d is  the  source  dimension.  The  corner  frequency,  fQ,  is 
defined  at  the  intersection  of  the  two  linear  segments 
which  represent  the  spectrum  and  is  shown  in  Figure  1. 

Following  Savage  (1966)  , whether  using  P-waves  or  S- 
waves,  the  displacement  signal  has  an  unidirectional  pulse 
character  and,  consequently,  the  displacement  spectrum  for 
any  of  the  body  phases  should  have  a maximum  at  zero  fre- 
quency and,  therefore,  the  displacement  spectrum  in  the 
frequency  range  between  o < f < fQ  should  have  zero  slope. 
This  behavior  of  the  spectrum  of  the  body  phases  in  the 


range  of  frequencies  less  than  the  peak  frequency  is  then 
conditioned  by  the  long-time  displacement  conduct  of  the 
source  and  is  very  probably  due  to  the  offset  on  both 
sides  of  the  faulting  surface.  The  expected  spectrum  is 
illustrated  in  Figure  1.  In  the  case  of  an  explosion/  it 
is  not  probable  that  permanent  differential  movement  takes 
place  beyond  the  equivalent  elastic  cavity,  or,  in  other 
words,  there  is  no  offset  as  in  t.he  case  of  an  earthquake 
with  a faulting  mechanism  and,  therefore,  the  slope  of 
the  displacement  spectrum  of  the  P or  S phases  will  be 
positive  for  frequencies  smaller  than  the  fQ.  This  is 
illustrated  in  Figure  2. 

For  frequencies  greater  than  f0,  the  spectral  dis- 
placement  will  reflect  the  rise  time  of  the  source  function 
and  is  expected  to  be  very  much  the  same  for  both  types  of 
seismic  events.  These  two  characteristics  of  the  amplitude 
spectra,  the  peak  frequency  and  the  slope  of  the  amplitude 
spectra  in  the  low  range  of  frequencies  (f  < fQ) , are 
potential  tools  in  the  discrimination  between  natural  and 
manmade  seismic  events. 

None  of  the  theoretical  amplitude  spectra  corres- 
ponding to  the  existing  models  for  an  explosive  source 
have  peaks  and  valleys  in  counterpoint  with  those  from 
actual  events,  a fact  that  indicates  the  existence  of 
modulations  (Knopoff  and  Pilant,  1963) . The  causes  of 
these  modulations  can  be  numerous  but  among  the  most 
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cffoctive  arc  the  difference  in  the  path  due  to  the  layered 
nature  of  the  earth's  crust,  and  the  reverberation  of  the 
near-station  and  the  near-source.  In  the  case  of  a mul- 
tiple explosion  ^ith  spatial  separation  of  the  shot  points 
and/or  firing  time  delay,  a modulation  effect  would  also 
result  from  the  geometry  or  time  delays. 

The  last  modulating  effect  of  those  mentioned  above  is 
the  one  that  has  received  primary  attention  in  this  work. 
The  multiplicity  of  a nuclear  explosion  will  be  theoreti- 
cally investigated  in  several  possible  situations  under 
the  following  five  assumptions : 

1.  The  maximum  order  of  multiplicity  is  two. 

2.  Both  explosions  are  considered  to  be  of  the  same 
order  of  magnitude. 

3.  The  earth  and  the  recording  stations  behave  as 
linear  functions  in  terms  of  the  radiated  energy. 

4.  The  focuses  of  all  nuclear  explosions  are  suf- 
ficiently close  that  the  near-source  conditions 
are  equivalent. 

5.  The  source  area  is  transparent  to  propagation 
following  detonation. 

Let  e(t)  and  s(t)  be  the  impulse  response  of  the 
earth  and  recording  station,  respectively,  and  r1(t),  the 
time  function  that  represents  the  recorded  seismogram. 


It  follows  that: 
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rl(t)  = P(t)  * c(t)  * s(t)  (2) 

whore  p(t)  is  the  mathematical  expression  in  the  time 
domain  of  the  signal  emitted  by  the  source,  and  * denotes 
convolution . 

For  a double  explosion,  p (t)  can  be  expressed  as 
follows : 

p(t)  = f(t)  + f(t-ATx)  - Kf(t-AT2)  - Kf(t-AT3)  (3) 

where  f(t)  is  the  mathematical  expression  of  the  signal 
emitted  by  a single  explosion;  f(t-AT1),  accounts  for  the 
signal  emitted  by  the  second  explosion  with  a delay  time 
Aj.^;  Kf(t-AT2)  an<^  Kf(t-AT3)  are  the  time  functions  that 
represent  the  reflected  pulse  or  waves  in  the  free  surface 
of  the  first  and  second  explosions,  respectively;  K is  the 
scaling  factor  and  AT2  and  AT3  are  the  delay  times  such 
that  AT3  = AT2  + AT^ . 

Therefore,  combining  (2)  and  (3),  it  follows: 

ri(t)  = [f  (t)  + f(t-AT2)  - Kf(t-AT2)  - Kf  (t-A'T3)  ] * 

e(t)  * s(t)  (4) 

and  equation  (4)  in  the  Fourier  transform  domain  is  given 

by: 

= F(u)  *E(w)  *S  (u>)  (1  + e3UiiTl  - Ke^U^T2)  - Ke^W^T3)  (5) 

where  E(w)  and  S(w)  are  the  transfer  functions  of  the 
earth  and  the  recording  station;  R1(w)  and  F(w)  are  the 
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Fotirier  transforms 
angular  frequency. 


of  r^(t)  and  f(t);  and  w is  the* 


The  amplitude  spectra  of  R^ui)  is  then: 

,K1(w)  J - | F (oj)  | • j E ( oi)  | • J S (to)  | [2  + 2K2  + (1  + K2  - 
2K  cos  (0AT2)  *2  cos  ojAT^  - 4K  cos  ioAT2]^2 


The  factor 


[2  •»  2K2  + (1  + K2  - 2K  cos  <oAT9)  2 cos  wAT  - 

1 

4K  cos  toAT-5  i1/2 

in  (6)  is  the  modulating  factor  of  the  amplitude  spectra 
|Rl(io)  | due  to  the  second  explosion  and  the  two  pp  waves. 
It  should  be  noted  that  the  absence  of  AT3  in  the  modu- 
lating function  is  due  to  the  linear  relation  that  has 
been  assumed  between  AT3,  AT2  and  AT-^. 

In  the  case  of  a single  explosion,  the  other  effects 
being  equal 


r2  ( t)  = If(t)  - Kf(t-AT2H  * e(t)  * s(t) 
and  tne  Fourier  transform  of  (8): 

R2(w)  = F(u)) -E(oj) -S(u>)  (1  - Ke3wAT2j 


+ K2  - 2K  cos  u)AT2]1/2 


and  the  amplitude  spectra  is: 
lR2<“H  = |P(w)  | • JE{«)  J • jS(«)  | [1 


(10) 
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From  (10) , it  can  be  seen  that  a modulating  effect  also 
exists  and  is  due  to  the  reflected  waves  pP  in  the  free 
surface . 

As  initially  pointed  out  in  this  section,  the  near- 
source and  the  near-station  reverberations  will  have 
modulation  effects  upon  the  amplitude  spectra  which  are 
rather  difficult  to  eliminate  using  the  record  of  only 
one  shot  (whether  multiple  or  single) . An  alternative 
procedure  to  overcome  this  difficulty  is  to  use  the  spec- 
tral ratio  of  two  explosions  whose  focuses  lie  at  the 
same  location  and  which  have  been  recorded  at  the  same 
stations  with  the  same  equipment.  Under  these  conditions, 
the  effect  of  the  near-source  and  the  near-stations  will 
be  the  same  on  both  amplitude  spectra  and  therefore 
cancelled  in  the  ratio.  This  method  has  been  used  by 
Healy  and  others  (1971)  and  King  and  others  (1972,  1974). 

In  the  case  of  two  double  explosions,  the  amplitude  spectra 
ratio  will  have  the  form: 

[2  + 2K2  + (1  + K2  - 2K  cos  u)AT0)2  cos  wATn  - 
RATIO  (w)  = = ± — 

[2  + 2K2  + (1  + K2  - 2K  COS  wAT^)  2 cos  wA^  - 
1/2 

4K  cos  o)AT23 

7-rjz  (ID 

4K  cos  wAT^]1'^ 

In  the  case  of  a multiple  detonation  and  a single  one,  the 
corresponding  ratio  will  be: 
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ratio  ,w)  = 


[2  + 2K2  + (1  + K:  - 2K  cos  wAT2) 2 cos  - 4K  cos  ^V2I 

(12) 


1/2 


1 + K2  - 2K  cos 


OiAT2) 1/<2 


When  both  explosions  arc  single,  the  corresponding  ratio  of 
the  amplitude  spectra  will  have  t..e  following  form: 


RATIO  (u) 


(1  + K2  - 2K  cos  uAT2)1//2 
(1  + K2  - 2K  COS  wAT^)1  2 


(13) 


Correlation 

The  cross-correlation  function  is  a coherency  measure 
that  gives  in  a quantitative  manner  the  similarity  between 
two  sequential  sets  of  data.  For  data  in  a continuous 
form,  the  cross-correlation  function  $XY  between  two  time 
series  X(t)  and  Y(t*  is  given  by: 

1 T 

$VV(T)  = T / X ( t)  Y ( t + T)UT  ^14) 

AX  0 

where  x is  the  time  lag.  When  the  data  are  in  discrete 
form,  the  cross-correlation  is  given  as. 


Z ^XNX-J  ” ^ ^YJ  ~ 
J=1 


V.Y  “ NX  2 NY  , _|  2 1 1/2 

1 z lxT  - x|2  2 i*T  - Yi  I 

J=1  J J=1  J 


(15) 


where:  NX  = the  number  of  points  in  data  set  X 
NY  = the  number  of  points  in  data  set  Y 
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The  cross-correlation  function  is  used  in  this  study 
as  a measure  of  similarity  between  two  seismic  signals 
recorded  at  the  same  station  from  events  which  are 
theoretically  located  at  the  same  point.  Due  to  the  dif- 
ference in  the  source  mechanism  of  shallow  earthquakes  and 
nuclear  explosions,  the  signal  shape  and  its  frequency 
content  produced  by  these  two  different  sources  will 
differ  somewhat.  It  then  appears  reasonable  to  expect  a 
higher  value  of  the  cross-correlation  for  the  zero  lag 
when  comparing  two  nuclear  shots  than  when  an  earthquake 
and  a nuclear  detonation  are  correlated.  By  the  same 
line  of  reasoning,  a higher  correlation  will  be  expected 
for  two  double  detonations  than  when  a double  and  a single 
one  are  compared. 

Maintaining  the  same  assumption  for  the  nuclear  explo- 
sions enumerated  on  the  amplitude  spectra  theory,  let  us 
examine  in  analytical  form  the  resulting  cross-correlation 
function  for  several  cases. 

Let  r(t)  be  the  time  function  that  represents  the  P 
coda  of  an  earthquake's  recorded  signal  at  some  station 
composed  of  the  direct  P-waves  r^t)  and  the  pP  phase 
kr^(t  - AT^) , such  that: 
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r(t)  - r1(t)  - kr^tj  - AT^)  (16) 

where  AT^  is  the  time  delay  of  the  pP  waves,  and  k,  a 
scaling  factor. 

In  a similar  manner,  if  q1(t)  and  kq^t  - At^) 
represent  the  direct  P and  the  pP  phases,  respectively, 
from  a single  nuclear  explosion,  its  recorded  signal 
qs(t)  can  be  represented  as  follows: 

c3s(t)  = qi<t)  - kq1(t  - ATX)  (17) 

and  assuming  that  AT2  is  the’  delay  time  for  a second 
explosion  with  characteristics  identical  to  the  first  one, 
the  time  expression  qd(t)  of  a recorded  signal  for  a 
double  shot  will  be  represented  as  follows: 

qd  = qi(t)  - ^(t  - ATx)  + qx(t  - AT2)  - kqx(t  - ATX  - AT2) 

(18) 

Now  consider  the  correlation  function  for  the  following 
cases : 

1.  Earthquake  and  single  nuclear  explosion 

2.  Earthquake  and  double  nuclear  explosion 

3.  Single  nuclear  explosion  and  double  nuclear 
explosion. 

Earthquake-single  shot: 

1 T 

<j>rqf((T)  = ^ / [rx(t)  - krL(t  - ATX ) ] 

[qx(t  + T)  - kq1(t  + T - AT^ldt 


(19) 


1 
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?rqg  = (1  + k2)*e(T)  - k[e(x  + AT^)  + e(x  - AT1)  ] (20) 

where 

1 T 

e(x)  = m / r,  (t)  *q.  (t  + x)dx 
0 1 

Earthquake-double  shot: 

1 T 

0rqd(t)  = - / [r^t)  - kr^t  - AT1)][q1(t  + x)  - 

kq1(t  + x - AT^)  + q1(t  + x - AT2)  - kq1(t  + T - A^  - AT2)]dx 

(21) 

2 

<?rqd  = (1  + k )e(x)  - k[e(x  + AT^)  + e(x  - AT^)  + 
e(x  + AT1  + AT2)  + e(x  + AT2  - AT^) ] + 

(1  + k2) e (x  + AT2)  (22) 

Single  shot-double  shot: 

1 T 

<*>qsqd(T)  = Tq^  ^(t)  " ^(t  - AT1)](q1(t  + T)  - 
kq-j  (t  + T - AT-^ ) + q1(t  + T - AT2)  - 

kqi(t  + X - AT-j_  - AT2)]dx  (23) 

$q  <3,(t)  = (1  + k2)a(x)  - k[a(x  + AT-.)  + a(x  - AT-,)  + 

S a J- 

a(x  + AT^.  + AT2)  + a(x  + AT2  - AT^  ] + (1  + k2)a(x  + AT2) 

(24) 


♦ 


a (t) 


= V 


qn  (t)q,  (t  + x)dx 


where 
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Thorc  arc  three  other  trivial  cases:  earthquake- 

earthquake,  single  shot-single  shot  and  double  shot- 
double  shot,  for  which  obviously,  due  to  our  theoretical 
frame,  the  cross-correlation  function  is  coincident  with 
the  auto-correlation  function,  and,  consequently,  the 

value  of  the  correlation  function  for  zero  lag  is  +1  for 
all  of  them. 

From  expressions  (20)  and  (22),  the  following  rela- 
tion can  be  written: 

^rqd  = ^rqs  “ + ATX  + AT2)  + c(t  + AT2  - A?x)  ] + 

(1  + k2 ) e ( t + AT2)  (25) 

which  clearly  indicates  the  difference  that  exists  in  the 
value  for  any  t of  the  cross-correlation  function  for  the 
earthquake-singie  shot  and  the  earthquake-double  shot 
cases. 

The  difference  of  the  cross-correlation  function, 
that  results  when  comparing  single-shot  and  double-shot 
with  the  first  two  cases  presented,  is  obvious  by  inspec- 
tion of  the  expressions  (20),  (22),  and  (24).  Expressions 
(20)  and  (22)  both  contain  e(x),  which  is  the  cross- 
correlation  of  the  direct  P-waves  time  function  from  the 
earthquake  and  the  direct  P-waves  time  function  of  the 
single  nuclear  shot;  and  (24)  contains  the  auto-correlation 
of  the  P-waves  time  function  of  a single  shot. 


Theoretical  Source  Models 


The  source  dimension  for  the  Massachusetts  Mountain 
earthquake  was  calculated  using  a relation  similar  to  that 
given  by  Brune  (19  70,  1971)  and  in  which  the  velocity  of 
the  P-waves  has  been  substituted  for  the  S-waves  velocity. 
The  actual  expression  used  was: 


r 


2.34v 

2uf^ 


(26) 


where  v is  the  velocity  of  the  P-waves  in  the  immediate 
vicinity  of  the  source  and  fQ  is  the  corner  frequency  of 
the  P-waves  spectrum.  This  approach  fcr  the  estimation 
of  the  source  dimension  has  been  utilized  by  Wyss  and 
others  (1971)  and  Hanks  and  Wyss  (1972).  While  widely 
used,  in  application  the  validity  of  (26)  is  dependent 
upon  the  following  assumptions : 

1.  The  earthquake  is  produced  by  the  dislocation  of 
a strike-slip  fault. 

2.  The  dislocation  is  produced  by  a tangential 
stress  applied  to  the  interior  of  the  fault  area, 
which  is  stimulated  in  an  infinite,  homogeneous, 
perfectly  elastic  and  isotropic  medium. 

3.  The  rupture  is  initiated  instantaneously  and 
simultaneously  over  the  total  fault  area. 

4.  The  rupture  velocity  is  always  subsonic. 
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5.  During  the  time  that  the  rupture  takes  place,  the 
fault  surface  reflects  perfectly  the  shear  waves, 
or  in  other  words,  each  block  of  the  fault 
behaves  as  a liquid  with  respect  to  the  other. 

6.  The  stress  pulse  only  produces  shear  waves  which 
propagate  in  a direction  perpendicular  to  the 
fault  surface. 

Sharpe's  (1942)  model  for  an  explosive  source  was 
used  for  the  estimation  of  the  source  dimension  of  the 
nuclear  detonations.  In  Sharpe's  model,  the  corner  fre- 
quency is  related  to  the  P-waves  velocity  v of  the  media 
near  the  source  and  the  source  radius  r as  follows: 

' fQ  = /2  v/37tr  (27) 

The  validity  of  this  equation  is  dependent  upon  the  fol- 
lowing assumptions: 

1.  The  explosion  takes  place  in  a spherical  cavity. 

2.  The  medium  surrounding  the  spherical  cavity  is 

, infinite,  homogeneous  and  perfectly  elastic. 

3.  No  net  displacement  occurs  beyond  the  equivalent 
elastic  cavity. 

4.  An  arbitrary  pressure  function  of  time  is 
applied  uniformly  to  the  interior  of  the  cavity. 
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Geological  Setting 

The  Nevada  Test  Site  (NTS)  lies  between  the  east 
boundary  of  the  Cordilleran  eugcosyncline  and  the  east 
boundary  of  the  miogeosyncline , close  to  the  Las  Vegas 
Valley  shear  zone  (Figures  3 and  4) , within  the  Basin 
and  Range  province. 

The  rocks  that  constitute  the  Late  Precambrian  and 
Early  Cambrian  are  typically  clastic,  chert  being  the 
most  representative  lithologic  unit.  Silurian  and  Lower 
Devonian  rocks  are  mainly  dolomites  with  depositioaal 
characteristics  of  broad  shallow-subtidal,  intertidal 
environments  as  documented  by  Matti  and  others  (1974) . 

The  Roberts  Mountain  thrust  plate,  formed  northwest  of 
the  NTS,  was  the  source  of  part  of  the  clastic  material 
that  constitutes  the  Eleana  Formation  deposited  during 
Late  Devonian  and  Mississippian.  The  rest  of  the  Paleo- 
zoic is  characterized  by  carbonate  deposition. 

Mesozoic  sediments  and  volcanic  rocks  are  almost 
absent  in  the  area  of  study,  being  represented  only  by  a 
few  granitic  plutons  dated  in  that  period. 

Ash-flow  tuffs,  mainly  of  rhyolitic  composition 
dominated  the  lithology  of  the  more  than  3,700  foot  sec- 
tion which  constitutes  the  Tertiary  and  has  been  used  as  a 
test  media  for  many  of  the  underground  nuclear  detonations 
carried  out  in  the  NTS.  Given  in  Table  1 are  the  princi- 
pal Tertiary  volcanic  units  at  NTS  and  the  surrounding 


L 


Figure  3.  Physiographic  provinces  of  the  Nevada 
California,  and  Utah  area. 
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Figure  4.  Index  map  of  the  geologic  and  physiographic 
features  of  the  Nevada  Test  Site  and  sur- 
rounding area.  (modified  from  Ekren,  1968, 

p.  12). 


♦ 
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TABLE  1 


Principal  Tertiary  Volcanic  Units  at  NTS,  K/Ar  Dates, 
Areal  Extent,  and  the  Estimated  Volume  of  Each  Unit 


Volume 

Lateral 

In 

Extent, 

Cubic 

Age 

K/Ar 

Date 

Unit 

In  Miles 

Miles 

7.5 

m.y . 

Thirsty  Canyon  tuff 

80 

50 

11.0 

m.y. 

Timber  Mountain  tuff 

100 

500 

>< 

& 

Paintbrush  tuff 

70 

120 

< 

H 

Eh 

« 

EH 

13.8 

m.y. 

Belted  Range,  tuff 

60 

70-100 

17.5 

m.y. 

Fraction  tuff 

120+ 

500 

White  Blotch  Spring  tuff 

120+ 

500 

26.5 

m.y. 

Monotony  Valley  tuff 

120+ 

1000+ 

(Modified  from  E.  B.  Ekren,  1968) 
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arGa'  the  K/Ar  dates,  the  areal  extent,  and  the  estimated 
volume  of  each  unit. 

The  structural  features  of  the  NTS  are  those  of  the 
Basin  and  Range  province,  characterized  by  a north-south 
trend  of  mountain  ranges  and  intermontane  basins  due  to  an 
intense  faulting.  Shown  in  Figure  4 are  the  major  thrust 
faults  and  other  physiographic  features  of  the  NTS. 

All  the  seismic  events  used  in  this  study  have  been 
fired  or  have  occurred  in  the  Yucca  Flat  area,  which  is 
in  an  intermontane  basin  mainly  covered  with  Tertiary 
(volcanic)  ar.d  Quaternary  (alluvium)  deposits,  and  which 
is  bordered  by  mountains  consisting  of  Precambrian  and 
Paleozoic  rocks  and  a few  Mesozoic  intrusives  in  the 
northern  region.  A general  geologic  map  of  the  Yucca  Flat 
area  is  presented  in  Figure  5. 

A sound  knowledge  of  the  geology  of  the  area,  not  only 
where  the  observed  seismic  events  take  place  but  also  along 
the  propagation  path  separating  the  source  and  the  record- 
ing station,  is  of  prime  importance.  Much  of  what  is  con- 
tained on  a seismogram  is  the  result  of  the  geological 
structure  and  major  lithologic  units  that  constitute  the 
path  traveled  by  the  seismic  energy.  It  is  not  difficult 
to  realize  the  complexity  of  these  effects  and  the 
arduousness  of  their  evaluation.  Much  of  the  work  along 
this  line  for  the  NTS  area  has  been  dene  by  the  United 
States  Geological  Survey  (USGS) . Yacoub  and  others 
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□ 


Figure  5.  Generalized  geologic  map  of  Yucca  Flats  area 
with  epicenter  of  Massachusetts  Mountain 
earthquake  and  area  in  which  nuclear  events 
were  detonated.  (modified  from  Hinrichs, 
1968,  p.  242). 
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(1968,  1970),  have  studied  the  effects  of  the  subsurface 
geology  on  the  ground  motion  due  to  underground  nuclear 
detonations  at  the  NTS  by  using  computer  ray  tracing 
techniques  on  two  complex  two-dimensional  geological 
models.  Since  the  two  mentioned  models  include  the  sites 
of  the  Tonopah  and  Nelson  stations  (which  are  two  of  the 
data  sources  used  in  this  thesis) , a brief  review  of  the 
results  from  the  models  concerning  the  seismic  stations 
is  presented. 

Yacoub  and  others  (1970)  applied  the  ray  tracing  com- 
puter program  to  the  two  geological  models  indicated  by 
lines  AB  and  BC , respectively,  in  Figure  6 with  a seismic 
source  for  both  models  located  at  point  B which  is  coinci- 
dent with  drill  hole  UE20f.  Shown  in  Figures  7 and  8 are 
the  geological  model  and  the  program  output  that  correspond 
to  model  AB.  In  Figure  8,  the  high  ray  density  that  can  be 
observed  between  distances  of  330,000  and  350,000  feet 
bracketing  the  Tonopah  station  site  and  its  immediate  vicinity 
has  been  suggested  to  be  a possible  partial  explanation  of 
the  large  ground  motion  recorded  at  this  station  for  many 
events  in  the  NTS.  The  high  concentration  of  rays  (Figure 
8)  seems  to  be  controlled  mainly  by  the  lithological  unit 
(Figure  7)  to  the  Conrad  and  Moho  discontinuity.  An 
important  observation  (pointed  out  by  Yacoub  and  others, 

1970)  is  that  if  the  source,  rather  than  being  located 
in  UE20f,  were  located  at  a point  between  Tonopah  and 


Figure  7.  Geological  model  of  Line  AB.  (from  Yaucoub  and  others,  1970 
p.  591)  . 
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Ui:20f,  the  station  would  probably  fall  in  a "shadow  zone" 
(low  density  of  rays)  similar  to  that  of  the  Ralston 
Valley  and  Cactus  Range  shown  in  Figure  8.  This  point 
is  particularly  important  with  respect  to  the  emplacement 
of  the  stations  in  an  unmanned  monitoring  system. 

The  model  for  line  BC  is  not  as  detailed  as  that  for 
line  AB.  Hence, the  corresponding  geological  model  and 
output  have  not  been  presented  here.  It  appears  that 
the  majority  of  the  emerging  rays  at  the  Nelson  station 
have  undergone  many  refractions  in  the  crust  prior  to 
reaching  the  area  close  to  the  station  which  is  underlain 
by  intrusive  and  extrusive  rocks  with  a small  number  of 
elastic  interfaces. 

In  Table  2,  the  average  values  are  listed  for  the 
longitudinal  velocity,  density  and  Poisson's  ratio  for 
the  major  lithologic  units  in  the  NTS.  These  values  are 
used  in  this  study  and  were  taken  from  Yacoub  and  others 
(1970)  . 
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TABLE  2 


Longitudinal  Velocity,  Density  and  Poisson's  Ratio 
of  the  Major  Lithologic  Units  in  the  NTS 


Longitudinal 
Velocity 
km/ sec 

Density 

gr/cc 

Poisson 1 s 
Ratio 

Rock  Types 

1.2 

1.80 

.405 

Alluvium;  vitric 
bedded  tuff 

3.4 

2.30 

.330 

Welded  tuff 

2.4 

2.00 

.270 

Zeolitic  bedded  tuff; 
calcareous  sedimen- 
tary rocks 

4.6 

— 

2.40 

.300 

Silicic  volcanic  rocks 

5.5 

2.65 

.250 

Silicic  plutonic  and 
hypoabyssal  rocks  at 
shallow  depths; 
quartzite  and  argil- 
lite 

6.0 

2.70 

.250 

1 

Limestone 

6.7 

2.75 

.250 

1 

i 

Dolomite 

(Modified 

from  N . K . 

Yacoub  and 

others,  1970) 
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Data 

Nino  seismic  events  from  the  NTS  comprise  the  data 
used  in  this  study.  The  Massachusetts  Mountain  earthquake 
that  occurred  August  5,  1971,  on  the  NTS  between  the  Yucca 
and  Frenchman  Flats  (Fischer  and  others,  1972;  Rohrer  and 
Springer,  1972)  is  the  only  natural  event.  The  rest  are 
nuclear  explosions.  The  nuclear  explosion  BLENTON/THISTLE 
is  known  to  be  a double  detonation  and  DIDO  QUEEN,  a 
single  one;  the  remaining  shots  from  A to  F are  of  unknown 
multiplicity.  The  arrivals  from  an  earthquake  with  an 
epicenter  at  teleseismic  distance  from  the  shot  point  of 
Event  F occurred  almost  simultaneously  with  the  firing  of 
this  nuclear  detonation  and  both  events  were  recorded 

together.  In  Table  3,  the  data  for  all  events  are  sum- 
marized. 

The  seismographic  coverage  of  the  area  is  comprised 
of  the  Lawrence  Livermore  Laboratories  (LLL)  stations, 

ELK,  KNB,  LAC  and  MNV  and  the  Sandia  Laboratories  stations, 
BMN,  DAC,  ELY,  LEE,  NEL  and  TPH.  The  approximate  loca- 
tions of  these  stations  are  shown  in  Figure  9,  and  their 
geodetic  coordinates  given  explicitly  in  Table  4.  The 
azimuth  and  distance  of  each  event  from  BLENTON,  which 
serves  as  somewhat  of  a reference  for  this  study,  are 
also  given  in  Table  4.  The  recordings  at  the  LLL  and 
Sandia  stations  were  obtained  with  Sprengnether  and  Benioff 
instruments,  respectively.  The  amplitude  of  the  velocity 


Summary  of  the  Data  of  the  Seismic  Events  Used  in  this  Study 
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Events  A through  F are  all  located  at  distances  less  than  20  km  from 
BLENTOM/THISTI.E. 
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Figure  9.  Location  of  Nevada  Test  Site  and  Lawrence 
Livermore  (ELK,  KNB , MNV,  and  LAC)  and 
Sandia  (BMN , DAC,  ELY,  LEE,  NEL,  and  TPH) 
stations . 
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response  for  each  of  these  instruments  is  given  in  Figure 
10.  A summary  of  the  records  used  in  this  study  is  given 
in  Table  5.  The  records  indicated  by  this  table  were 
selected  on  the  basis  of  record  quality  and  availability. 


L 


Simmary  of  the  Records  Used  in  this  Study 


-38- 


£ 

u;  U4 

w 


x 


x 


2 

w ca  , x 

« ! 


xxxxxxxx 


2 < 1 
> I 

.-M  , 


X X X X X 


Gt 


b* 

•73 


c-* 

\ 


-4 

A 


■P 

O 

03 

2 

o 

<3 

(.3 


o 

•5  n 

-p  rr 

r* 

3 D 
o u 


5 w 


-3 

W 


m 


u 

5 


> 

g 


fli 


0 

>« 

W 

22 

< 

u: 

w 

e* 

Q 

U 

*1 

z 

H 

J 


1 


-39- 


Digital  Analysis 

Approximately  30  seconds  of  the  P-wave  and  P-wave 


coda  for  all  records  were  digitized  with  a Computer  Equip- 
ment Corporation  digitizer.  The  initial  digitization  was 
performed  at  unequally  spaced  intervals  in  order  to  obtain 


maximum  information  from  the  records  while  minimizing 
digitizing  efforts.  The  unequally  spaced  data  was  then 
lintarly  interpolated  at  equal  0.05  second  increments  to 
facilitate  computer  processing.  This  uniform  rate  of 
sampling  yields  a folding  (or  Nyquist)  frequency  f^  of 


f - 1 1 

rN  “ 2At  “ 2 ( .05) 


10  Hz. 


(28) 


On  the  basis  of  electrical  and  material  considerations, 
this  Nyquist  frequency  is  satisfactory. 

The  amplitude  spectra,  amplitude  spectral  ratios  and 
the  natural  logarithm  of  the  amplitude  spectra  of  the  data 
were  calculated  using  a computer  program  which  performs  the 
operations  schematized  in  Figure  11.  The  tapering  opera- 
tion was  applied  to  smooth  the  ends  of  the  data  vector  and 
is  equivalent  to  the  application  of  a tapered  time  window 
to  the  data  set.  The  taper  used  is  given  by 


i(l  + cos  2 fry  2 ) , J * 1,3 


T(J) ■-  < 


, J - 4,  NP-3 


(29) 


+ eosliN-L..^.  3 " J) , J - NP-2,  NP 
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where  NP  is  the  number  of  points  in  the  input  data  set. 

The  Fourier  transformation  was  accomplished  with  a FFT 
(Fast  Fourier  Transform)  package  due  to  Cooley  and  Tukey 
(1965)  . 

The  cross-correlations  were  calculated  in  the  time 
domain  with  a computer  program  aue  to  Taylor  (1974) . Thr 
cross-correlation  program  is  conventional  in  all  respects 
and  utilizes  the  following  definition  of  the  cross— correla- 
tion at  lag  NAt. 


NX-N  _ ' 

{Jl,  XNX-N-K  ^ X yK  “ Y) 

$XY(NDT)  = 

NX  _ , NY  , w_ 

I 2 [X  - X]2  Z [V  - Y)2\1/2 
K=1  K K=1  K 

where 

^XY  = cross“correlation 
NX  = the  number  of  points  in  lata  set  X 

NY  = the  number  of  points  in  data  set  Y 

X = the  mean  of  data  set  X 

Y = the  mean  of  data  set  Y 


(30) 
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Results 

Spectral  Analyses 

For  an  = 5 explosion  with  an  assumed  source  radius 
of  0.6  km,  detonated  in  material  with  a P-wave  velocity  of 
2.5  km/sec,  equation  (27)  yields  a peak  frequency  of  0.63 
Hz.  The  above  parameters  are  typical  of  the  events 
employed  in  this  study  and  peak  frequencies  on  the  order 
of  0.63  Hz  are,  therefore,  expected.  From  Figure  10  it  is 
evident  that  this  frequency  falls  within  the  pass-band  of 
the  Sprengnether  system,  used  by  the  LLL  stations,  and  no 
instrument  correction  is  required.  For  the  Benioff  sys- 
tem, used  by  the  Sandia  Laboratories  stations.  Figure  10 
indicates  the  expected  frequency  of  approximately  0.63  Hz 
falls  outside  the  -6  db  point  of  the  pass  band.  The  slope 
of  the  amplitude  spectra  at  frequencies  less  than  fQ  for 
the  Sandia  stations  would  clearly  be  determined  by  the  low 
frequency  roll-off  of  the  sensing  system  and  would  be 
unrelated  to  signal  characteristics.  The  system  roll-off, 
thus,  invalidates  the  use  of  low  frequency  displacement 
spectra  slope,  for  the  Benioff  system,  as  a criterion  of 
differentiation.  While  corrections  for  instrument  response 
in  these  cases  may  have  been  possible,  it  was  felt  that 
the  corrections  over  the  required  part  of  the  band  would 
enhance  noise  to  an  unacceptable  level.  For  this  reason 
only  the  LLL  stations  with  their  broad  pass-bands  were 
employed  in  the  spectral  characteristics  portion  of  this 
study . 
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In  order  to  clearly  define  the  corner  frequencies  in 
the  displacement  spectra,  broad  band  frequency  information 
is  required.  The  amplitude  spectra  of  all  events  con- 
sidered here  were  determined  for  frequencies  from  DC  to 
10  Hz.  This  band,  however,  exceeds  the  signal  band  and 
the  spectra  are,  therefore,  presented  over  the  range  from 
.04  Hz  to  3 Hz. 

The  natural  logarithm  of  the  amplitude  spectra, 
uncorrected  for  instrument  effects,  versus  the  natural 
logarithm  of  the  frequency  for  five  typical  events  are 
shown  in  Figures  12  to  16.  The  spectra  for  all  other 
remaining  event-station  combinations  are  included  in 
Appendix  A.  The  single  nuclear  event,  DIDO  QUEEN,  was  not 
used  here  since  LLL  records  were  not  available  for  this 
event.  The  corner  frequencies  obtained  from  these  spectra 
are  summarized  in  Table  6.  The  apparent  source  radii  were 
calculated  from  the  observed  corner  frequencies  for  spectra 
and  assumed  P-wave  velocities  of  3.4  km/sec  for  the  Massa- 
chusetts Mountain  earthquake,  1.4  km/sec  for  the  nuclear 
explosions  fired  in  alluvium,  and  2.4  km/sec  for  those 
detonated  in  tuff.  The  calculated  source  radii  for  all 
events  are  summarized  in  Table  7.  The  average  radius  for 
each  event,  obtained  from  the  data  of  Table  7,  is  shown 
with  the  corresponding  depth  and  magnitude  of  the  event 
in  Table  8.  It  is  evident  from  Table  8 that  although  the 
magnitude  of  the  Massachusetts  Mountain  earthquake  is  less 
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Figure  12.  Displacement  spectrum  of 
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TABLE  8 

Average  Source  Radius,  Depth,  and  Magnitude  (mb) 
of  Massachusetts  Mountain  Earthnuake  and  Nuclear  Explosions 
BLEN TON /THISTLE , A,  B,  C,  D,  E and  F 


Event 

Average  Source 
Radius  (km.) 

Depth 
(km. ) 

Magnitude 

(mb) 

Massachusetts  Mountain 

1.1 

4.6 

4.3 

BLENTON/THISTLE 

.47 

.56 

5.3 

EVENT  A 

.28 

.24 

4.7 

EVENT  3 

.50 

.56 

5.6 

EVENT  C 

.36 

.31 

5.0 

EVENT  D 

.47 

.44 

5.3 

EVENT  E 

.23 

.39 

5.3 

F (as  an  explosion) 

1.00 

.24 

r> 

• 

F (as  an  earthquake) 

3.33 

o 

• 

•> 

• 

-52- 


than  the  magnitude  of  the  nuclear  explosions,  its  average 
source  radius  is  clearly  greater.  It  was  previouslv  noted 
that  Event  F consisted  of  a natural  earthquake  arrival 
preceding  the  shot  arrival.  This  is  clearly  suggested  in 
Table  8,  indicating  a source  radius  too  large  for  a 
nuclear  event  of  this  magnitude . 

The  source  radii  versus  depth  and  magnitude  are 
plotted  in  Figures  17  and  18.  It  is  interesting  to  note 
that  for  all  explosions  with  the  exception  of  Event  E,  the 
average  source  radius  and  the  corresponding  depth  are  of 
the  same  magnitude.  The  same  relation  has  been  found  by 
Wyss  and  others  (1971)  for  the  MILROW  nuclear  explosion 
which  was  fired  at  a depth  of  1.2  km  in  a calculated  source 
radius  of  1.25  km.  These  authors  have  suggested  that  the 
coincidence  between  the  elastic  source  radius  and  the 
depth  of  the  event  could  be  due  to  a restraining  effect 
that  the  free  surface  exerted  upon  the  elastic  source 
dimension.  It  is  also  apparent  by  inspection  of  Figure  18, 
in  which  the  calculated  elastic  source  radius  of  events 
BLENTON/THISTLE,  A,  B,  C,  D and  E are  plotted  versus  the 
corresponding  magnitudes,  that  with  the  exception  of  Event 
E,  an  increase  in  the  magnitude  of  the  event  is  associated 
with  an  increase  of  the  calculated  equivalent  elastic 
source  radius. 

A difference  in  the  slope  of  the  displacement  spectra 
for  the  band  from  .20  Hz  to  the  corresponding  peak 
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Figure  17.  A/erage  source  radius  of  BLENTON/THISTLE , 
A,  B,  C,  D and  E events  versus  depth  of 
each  event. 
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Figure  18.  Average  source  radius  of  BLENTON/THISTLE, 
A,  B,  C,  D and  E events  versus  magnitude 
of  each  event. 
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frequency  for  the  Massachusetts  Mountain  earthquake  and 
the  nuclear  explosions  is  evident  from  Figures  12  to  16. 

An  inspection  of  the  remaining  displacement  spectra  in 
Appendix  A clearly  establishes  the  consistency  of  this 
difference.  The  mean  slope  of  the  earthquake  in  the  .04 
to  peak  frequency  band  is  near  zero,  while  the  slopes  of 
the  nuclear  shots  are  clearly  positive.  Event  F is  again 
an  exception  to  this  but,  as  previously  noted,  the  initial 
arrivals  for  the  event  were  the  result  of  an  earthquake. 

difference  in  the  oehavaor  of  the  displacement 
spectra  ror  long  periods  has  been  observed  by  Wyss  and 
others  (1971) . Those  authors  suggested  that  the  flatness 
of  the  displacement  spectrum  of  an  earthquake  at  low  fre- 
quencies is  indicative  of  a step  function  at  the  source, 
while  the  faster  deterioration  of  the  explosion's  dis- 
placement spectra  could  be  related  to  the  depth  of  the 
source  and  its  time  function.  On  the  basis  of  these 
limited  cases,  the  slope  of  the  low  frequency  part  of  the 
displacement  spectra  is  clearly  a differentiating  quantity. 

Special  attention  has  been  given  to  Event  F due  to 
the  ..act  that  it  is  a "mixture"  of  an  earthquake  that 
occurred  at  teleseismic  distance  and  a nuclear  detonation. 
It  is  evident  from  the  displacement  spectra  of  this  event 
(Figure  16)  that  its  low  frequency  behavior  is  that  of  a 
typical  earthquake.  The  source  dimension  for  Event  F was 
calculated  considering  its  dual  nature  as  a nuclear 
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explosion  and  earthquake.  The  results  are  given  in  Tables 
7 and  8.  Since  the  earthquake  occurred  at  teleseismic 
distance,  correction  for  attenuation  should  be  applied, 
and  the  3.33  km  obtained  for  its  source  radius  is  of  an 
approximate  nature.  Nevertheless,  its  calculated  source 
radius  clearly  indicates  an  earthquake.  The  source  dimen- 
sion obtained  by  considering  the  event  as  being  only  a 
nuclear  explosion  is  1.1  km  and  although  its  magnitude  was 
not  known,  this  radius  appears  to  be  too  large  for  an 
explosion  with  a depth  of  .24  km,  therefore  suggesting  an 
earthquake . 

With  the  techniques  considered  here,  it  does  not  appear 
possible  to  detect  the  presence  of  the  shot,  although  it 
should  be  noted  that  no  attempt  was  made  in  this  study  to 
extract  the  nuclear  explosion  from  the  record.  A possible 
explanation  of  the  earthquake-type  behavior  of  the  dis- 
placement spectra  of  Event  F can  be  made  by  assuming  the 
low  frequencies  result  primarily  from  the  earthquake.  In 
this  sense  the  event  should  be  considered  an  earthquake 
rather  than  an  underground  explosion.  On  this  basis,  the 
event  lends  significance  to  the  use  of  spectral  slope  and 
source  radii  as  methods  of  identification. 

A careful  examination  of  the  spectra  for  all  nuclear 
events  will  show  that  a strong  peak  occurs  consistently  at 
about  .16  Hz.  In  the  few  exceptions  which  do  exist,  a 
strong  peak  does  appear  in  the  neighborhood  of  .16  Hz. 
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This  peak  is  probably  unrelated  to  any  characteristic  of 
the  events  and  is  due  instead  to  the  seismic  background 
noise.  . The  existence  of  a peak  in  the  spectrum  of  ground 
noise  at  a frequency  of  .16  Hz  has  been  established  by 
Brune  and  Oliver  (1359)  and  Brune  (1971). 

Records  of  the  single  explosion  DIDO  QUEEN  from  the 
LLL  stations  were  unfortunately  not  available  and  compara- 
tive studies  of  the  behavior  in  terms  of  the  logarithm  of 
the  spectra  of  single  and  double  explosive  sources  could 
not  be  accomplished.  For  the  sake  of  completeness,  the 
logarithm  of  the  spectra  of  some  nuclear  detonations  from 
the  Sandia  Laboratories  stations,  BMN,  ELY,  NEL  and  TPH 
are  shown  in  Figures  19  to  22.  These  spectra  exhibit  a 
positive  slope  in  the  long  period  range,  but  due  to  the 
narrow  band  of  the  recording  equipment,  much  of  the  behavior 
of  the  spectra  in  the  lower  range  of  frequencies  is  condi- 
tioned by  the  amplitude  response  of  the  seismic  system. 
Therefore,  no  consideration  has  been  given  to  these  spectra 
in  the  calculation  of  the  source  radius. 

In  an  attempt  to  investigate  the  possibility  of  a 
multiple  explosive  nature  for  Events  A to  E,  the  expres- 
sions of  the  modulation  factor  given  in  the  Theory  section 
by  equations  (7)  and  (10)  were  calculated  for  various 
assumed  values  of  the  delay  times  AT^,  AT2,  and  AT3 . 

Typical  forms  of  the  resulting  modulation  function  are 
shown  in  Figure  23.  It  is  evident  from  this  figure  that 
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Figure  23.  Typical  form  of  the  theoretical  modulation 
functions.  (a)  single  explosion;  (b) 
double  explosion. 


the  modulation  function  is  highly  dependent  on  the  delay 
times  and  is  not  a slowly  varying  function.  The  theoreti- 
cal modulation  function  was  then  compared  to  the  actual 
amplitude  spectra  of  each  event.  Attempts  were  then  made 
to  match  the  peaks  and  troughs  of  the  modulation  function 
with  the  envelope  of  the  signal  spectra  by  trial  and  error 
changes  in  the  assumed  values  of  AT1#  AT2  and  AT^ . 

Examples  of  the  amplitude  spectra  for  which  a seemingly 
suitable  modulation  function  was  found  are  shown  in  Figures 
24  to  33.  In  the  remaining  cases  no  reasonable  modulation 
function  could  be  found.  In  Table  9,  the  times  of  the 
different  arrivals  are  given  and  also  their  character  found 
with  the  modulation  functions  which,  in  general,  appear  to 
be  consistent  with  those  found  by  Cepstrum  Analysis 
(Zimdars,  1974).  It  is  obvious  from  Figures  24  to  33  that 
complex  modulation  effects  are  contained  in  the  amplitude 
spectra.  As  discussed  in  the  Theory  section  of  this 
investigation,  the  causes  of  the  modulations  are  varied 
and  difficult  to  evaluate,  if  not  impossible.  This, 
coupled  with  the  first  order  approach  that  constitutes  the 
emi — rical  models  for  the  modulation  functions  used  in  this 
study,  makes  the  assumed  times  and  their  polarity  highly 
speculative.  Nevertheless,  the  results  of  this  approach 
complement  the  other  results. 

The  amplitude  spectra  considered  here  result  from  the 
combined  actions  of  the  source  function  and  propagation 
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Figure  27.  Amplitude  spectrum  of  Event  B (solid  line) 
and  modulation  function  (dotted  line)  for 
ELK  station. 
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Figure  30.  Amplitude  spectrum  of  Event  C (solid  line) 
nd  modulation  function  (dotted  line)  for 
KNB  station. 
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Figure  31.  Amplitude  spectrum  of  Event  D (solid  line) 
and  modulation  function  (dotted  line)  for 
ELK  station. 
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path  effects.  It  is,  therefore,  possible  that  spectral 
characteristics  attributed  to  the  source  actually  result 
from  propagation  path  effects.  In  the  Theory  section  it 
was  shown  that  forming  the  spectral  ratio  for  two  events 
would  remove  propagation  effects  and,  thus,  yield  addi- 
tional information  concerning  the  source  functions.  Typi- 
cal examples  of  the  calculated  spectral  ratios  are  shown 
in  Figures  34  to  36.  It  is  evident  from  Figure  34  that 
the  spectral  contents  of  the  signals  from  events  BLENTON 
and  D as  seen  at  Kanab  were  nearly  identical.  This 
would  suggest  an  equivalent  nature  for  the  two  events. 

The  spectral  ratio  for  the  same  two  events  as  seen  at 
Nelson  is  shown  in  Figure  35.  On  the  basis  of  this  figure 
it  appears  that  the  spectral  contents  of  the  two  events 
are  unrelated.  In  general,  the  spectral  ratios  were 
characterized  by  this  lack  of  consistent  results  and  the 
strong  dependence  on  frequency  which  is  obvious  in  Figure 
36. 

Using  equations  (11)  and  (12)  , attempts  were  made  to 
determine  the  theoretical  envelopes  of  calculated  spectral 
ratios.  Because  of  the  number  of  parameters  involved,  no 
reasonable  results  were  obtained. 
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Figure  34.  Amplitude  spectral  ratio  of  event 
BLENTON/THISTLE  to  Event  D for  KNB 
station. 
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Correlation  Analyses 

In  keeping  with  the  theoretical  development  previously 
considered , it  is  desirable  to  consider  the  correlation  of 
the  following  four  cases;  cross-correlation  of  natural 
events  with  nuclear  events,  known  single  explosion  with 
known  double  explosion,  shots  of  unknown  multiplicity  with 
single  detonation,  and  known  double  explosion  with  detona- 
tions of  unknown  multiplicity.  These  four  cases  are  repre- 
sented respectively  by  the  cross-correlation  of  all  nuclear 
explosions  with  the  Massachusetts  Mountain  earthquake;  the 
cross-correlation  of  event  B LENTON /THISTLE  with  single 
explosion  DIDO  QUEEN;  the  cross-correlation  of  DIDO  QUEEN 
with  explosions  B,  D and  E;  and  finally  the  cross-correla- 
tion of  BLENTON/THISTLE  with  events  A to  E. 

The  actual  cross-correlations  are  time  dependent 
functions  having  durations  of  approximately  60  seconds.  A 
typical  example  is  shown  in  Figure  37  for  20  percent  of  the 
lags  or  +6  seconds.  In  view  of  the  geographical  proximity 
of  all  the  events,  however,  the  significant  part  of  the 
correlation  should  exist  near  a lag  of  zero  seconds.  On 
this  basis  the  correlation  functions  presented  here  are 
given  for  only  7 percent  of  the  total  lags  or  +2  seconds. 
Since  the  correlations  are  conducted  in  the  time  domain 
and  the  two  events  are  recorded  at  the  same  station,  the 
narrow  band-pass  of  the  Benioff  seismic  system  used  in  the 


-80- 


Sandia  Laboratories  stations,  has  no  detrimental  conse- 
quences and,  therefore,  the  records  of  these  stations  may 
be  used. 

The  results  of  the  correlation  of  the  explosions  with 
the  Massachusetts  Mountain  earthquake  are  presented  in 
Figures  37  to  43.  By  inspection  of  these  figures,  one  can 
notice  that  in  general  an  increase  of  the  correlation 
occurs  for  values  of  the  time  lag  other  than  zero.  It  is 
also  interesting  to  note  the  poor  correlation  obtained  for 
Event  F despite  its  partial  earthquake  nature.  This  prob- 
ably results  from  the  difference  in  distance  from  the 
respective  epicenters  to  the  recording  stations  and  also 
is  due  to  the  fact  that  Event  F is  not  a "pure"  earthquake 
but  a "mixture"  of  an  earthquake  and  a nuclear  detonation. 
The  correlation  values  at  zero  lag  for  the  Massachusetts 
Mountain  earthquake  and  each  of  the  nuclear  explosions  are 
summarized  in  Table  10.  From  this  table  it  can  be  seen 
that  in  general  a wide  variation  exists  in  the  values  of 
the  cross-correlation  at  zero  lag.  Despite  this  wide 
variation,  the  maximum  value  obtained  was  0.2'J  and  typical 
values  were  on  the  order  of  0.1. 

The  cross-correlation  function  of  event  BLENTON/ 
THISTLE  with  the  single  explosion  DIDO  QUEEN  is  illustrated 
in  Figure  44,  from  which  it  is  apparent  that  a low  correla- 
tion exists.  The  correlation  value  at  zero  lag  of  .OS, 
given  in  Table  11 , is  much  i^wer  than  the  expected  value 


Figure  38.  <~ross-correlations  of  Massachusetts  Mountain  earthquake  with 
Event  B for  stations  (a)  ELK,  (b)  LAC,  and  (c)  MNV. 

Abscissa:  Time  lag  in  sec.  Ordinate:  Correlation 


Figure  40.  Cross-correlation  of  Ma 
Event  D for  stations  (a 
Abscissa:  Time  lag  in 


^husetts  Mountain  earthquake  with 
JB , (b)  LAC,  and  (c)  MTJV. 


Ordinate:  Correlation 
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TAD  LE  10 

Correlation  Values  at  Zero  Lag  of  Massachusetts  Mountain 
Earthquake  with  Nuclear  Explosions  BLENTON/THISTLE , 

A,  B,  C,  D,  E and  F for  the  LLL  Stations 


ELK 

KNB 

LAC 

KNV 

M& ss . Mt • 
BLEXTOX/THISTLE 

.09 
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Mass.  Mt. 
Event  A 
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.12 
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o 

• 

i 

.24 

Mass.  Mt. 
Event  C 

.29 

.21 

.05 

Mass.  Mt. 
Event  D 

-.12 

.10 

-.17 

.13 

Mass.  Mt. 
Event  E 

• 

o 

CO 

-.14 

.04 

'.ass.  Mt. 


Event  F 
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TABLE  11 

Correlation  Value  at  Zero  Lag  of  Single  Nuclear  Explosion 
DIDO  QUEEN  with  Double  Explosion  BLENTON/TIIISTLE  and 
Unknown  Explosions  B,  D and  E for  the  LLL  Stations 


Events 

Station 

Correlation 

DIDO  QUEEN  - BLENTON/THISTLE 

LEE 

• 

O 

CO 

DIDO  QUEEN  - EVENT  B 

DAC 

.23 

DIDO  QUEEN  - EVENT  D 

DAC 

.23 

d::do  queen  - event  e 


DAC 
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for  the  correlation  of  two  nuclear  detonations.  However, 
little  significance  can  be  given  to  the  results  of  only 
one  station.  The  results  of  the  correlation  of  Events 
3,  D and  E with  DIDO  QUEEN  are  also  shown  in  Figure  44. 

By  simple  inspection  of  this  figure,  it  is  apparent  that 
higher  correlations  exist  for  these  events  than  for 
BLENTON/TKISTLE . The  corresponding  correlation  values  at 
zero  lag,  which  are  summarized  in  Table  11,  indicate  that 
higher  correlation  values  have  been  obtained  than  when  the 
same  events  B,  D and  E were  correlated  with  the  Massachu- 
setts Mountain  earthquake.  Nevertheless,  the  same  assertion 
that  was  made  regarding  the  correlation  of  BLENTON/THISTLE 
and  DIDO  QUEEN  applies  here  and  not  much  meaning  can  be 
derived  from  the  results  of  only  one  station. 

The  correlations  of  Events  A to  E with  the  double 
detonation  BLENTON/THISTLE  are  shown  in  Figures  45  to  52, 
ana  the  corresponding  zero  lag  values  are  itemized  in 
Table  12.  Of  particular  interest  are  the  correlation 
values  obtained  for  BLENTON/THISTLE  and  Event  B given  in 
Table  12.  These  values  are  on  the  order  of  0.6  and  indi- 
cate a strong  similarity  of  both  events.  The  values 
obtained  in  the  cross-correlation  of  BLENTON/THISTLE  and 
Events  D and  E for  the  stations  LAC  and  KNB  are  also  high. 
The  similarity  of  the  events  indicated  by  the  cross- 
correlation values  at  zero  lag  for  the  mentioned  cases  is 
r_rongly  supported  by  visual  comparison  of  their  records 


Figure  46.  Cross-correlation  of  Event  B with  the  double  detonation 
BLENTON /THISTLE  for  stations  (a)  TPH,  (b)  DAC,  (c)  LAC, 
and  (d)  NEL. 

Abscissa:  Time  lag  in  sec.  Ordinate:  Correlation 


double  detonation  BLENTON/THISTLE  for 
station  BHN. 


Abscissa:  Time  lag  in  sec. 

Ordinate:  Correlation 
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Figure  50.  Cross-correlatio  \ of  Event  D with  the 
double  detonation  BLENTON/THISTLE  for 
station  ELY. 

Abscissa:  Time  lag  in  sec. 

Ordinate:  Correlation 
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Figure  51.  Cross-correlation  of  Event  E with  the  double  detonation 
BLENTON/THISTLE  for  stations  (a)  KNB,  (b)  NEL,  (c)  DAC, 
(d)  ELY,  (e)  LAC,  and  (f)  TPH. 

Abscissa:  Time  lag  in  sec.  Ordinate:  Correlation 
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Fiqure  52.  Cross-correlation  of  Event  E with  tne 
double  detonation  BLENTON/THISTLE  for 
station  LEE. 

Abscissa:  Tine  lag  in  sec. 

Ordinate:  Correlation 
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TABLE  12 

Correlation  Values  at  Zero  Lag  of  Double  Explosion 
BLENTON/THISTLE  with  Unknown  Explosions  A,  B,  C,  D,  and  E 
for  LLL  and  Sandia  Laboratories  Stations 


BMN 

DAC 

ELY 

KNB 

LAC 
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TPH 
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.25 
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.57 
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.63 
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EVENT  C 
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.17 

.26 

.21 

.34 

.23 

BLENTON/THISTLE 
EVENT  D 

.23 

.47 

.55 

.33 

.25 

.42 

BLENTON/THISTLE 
EVENT  E 
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• 

.26 

.55 

.50 

.29 

.36 

.17 
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which  arc  shown  in  Figures  53  to  56.  The  correlation 
values  obtained  for  Events  A and  C for  all  the  stations 
are  lower  than  those  obtained  for  Events  R,  D and  E. 

These  results,  considered  in  the  theoretical  frame  pre- 
sented in  the  -Theory  -section  of  this  study,  suggest  a 
multiple  source  for  Event  B and  a single  explosion  for 
Events  A and  C.  No  indication  of  the  nature  of  Events 
B and  E is  provided  by  the  correlation. 

In  an  attempt  to  explain  the  variation  of  the  cor- 
relation values  obtained  for  these  events  at  different 
stations,  the  values  of  the  correlation  at  zero  lag  for 
each  event  have  been  plotted  versus  the  azimuth  of  the 
corresponding  stations.  The  correlation  values  for  the 
correlation  of  Events  A through  E with  3LENT0N/THISTLE 
at  each  station  are  shown  in  Figure  57.  On  the  basis  of 
this  figure,  it  is  apparent  that  for  the  majority  of  the 
events  considered  the  correlation  is  somewhat  dependent 
upon  the  azimuth.  In  general,  the  higher  values  appear 
in  the  neighborhood  of  200  degrees. 

In  Figure  58,  the  correlation  values  of  each  event  at 
each  station  are  plotted  versus  the  difference  in  depth. 
The  depth  difference  was  defined  as  the  difference  between 
the  depth  of  3LENT0N/THISTLE  and  the  depth  of  the  event 
being  considered.  It  can  be  seen  by  inspection  of  this 


figure  that,  in  general,  as  the  difference  in  depth  between 
the  sources  of  the  events  being  correlated  increases,  the 
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blenton/thistle  DAC  - VWB 


EYENT-B  DAC  - VWB 


Figure  53.  Analog  records  of  BLENTON/THISTLE  and 
B events  recorded  at  DAC  station. 


I 


BLENTON/THISTLE  TPH  - RWB 


EVENT-B  TPH  -RWB 


54.  Analog  records  of  BLENTON/THISTLE  and 
B events  recorded  at  TPH  station. 


Figure 


benton/thistle  lac  - mg 


EVENT-B  LAC  - MG 


EVENT-D  LAC  - MG 


EVENT-E  LAC  - MG 


Analog  records  of  BLENTON/THISTLE,  B,  D,  and 
E events  recorded  at  LAC  station. 
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Figure  56.  Analog  records  of  BLEN TON/THISTLE,  D,  and  E 
events  recorded  at  KN3  station. 
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Figure  58.  Correlation  values  at  zero  lag  of  A,  B,  C,  D and  E events  with 
BLENTON/THI STLE  explosion  versus  the  differences  in  depth 
between  each  event  and  BLENTON/THI STLE . 
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correlation  at  zero  lag  decrease a . This  depth  differc 
dependence  appears  to  limit  the  utility  of  correlation 
methods  for  identification. 


i ce 


Considering  the  positive  results  obtained  with  the 
limited  number  of  events  used  in  this  study  regarding  the 
different  behavior  of  the  slope  of  the  displacement  spectra 
in  the  lower  range  of  frequencies,  in  addition  to  the  pos- 
sibility of  this  test  being  performed  by  analogue  methods, 
further  research  should  be  conducted  to  provide  a more 
proper  evaluation  of  the  efficacity  of  this  test  and  give 
it  statistical  support.  Any  future  application  of  the 
slope  test  should  be  conducted  with  data  obtained  with  a 
broad  band-pass  seismic  system. 

Further  research  should  be  done  toward  the  investiga- 
tion of  particular  arrangements  of  multiple  explosions 
having  an  additive  effect  on  the  lower  frequencies  which 
might  cause  the  failure  of  the  slope  test. 

Based  on  the  differences  of  the  calculated  source 
radii  for  the  different  events  utilized  in  this  study, 
further  research  in  this  direction  is  recommended.  Although 
based  on  a limited  number  of  data,  the  positive  results  of 
the  correlations  used  in  this  study  suggest  that  these 
techniques  are  a potential  method  of  discrimination  between 
natural  and  manmade  seismic  events  and  identification  of 
multiple  nuclear  events.  Additional  investigation  on  a 
more  complete  set  of  data  than  that  used  in  this  research 
is  suggested  in  order  to  fully  evaluate  the  differentiating 
abilities  of  the  method. 
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APPENDIX  A 


DISPLACEMENT  SPECTRA  OF  MASSACHUSETTS  MOUNTAIN 


EARTHQUAKE  AND  NUCLEAR  EXPLOSIONS 


Figure  A1 . Displacement  spectrum  of  the  Massachusetts  Mountain  earthquake  for 
ELK  station  (loq  2.7,  Hz  = 0). 


spectrum  of  the  Massachusetts  Mountain  earthquake  for 
(logQ  2. 7,  Hz  = 0)  . 


LOG  FREQUENCY,  Hz. 

Figure  A4 . Displacement  spectrum  of  BLENTON/ THISTLE  event  for  LAC  station 
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Figure  A5 . Displacement  spectrum  of  A event  for  LAC  station 
(logG  2.7,  Hz  = 0) . 


Displacement  spectrum  of  B event  for  MNV  station 
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Figure  A8 . Displacement  spectrum  of  C event  for  KNB  station 
(loge  2.7,  Hz  = 0) . 
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Figure  All.  Displacement  spectrum  of  D event  for  ELK  station 
(log  2.7,  Hz  = 0) . 
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Figure  A13.  Displacement  spectrum  of  D event  for  LAC  station 
(logp  2.7,  Hz  = 0) . 


Figure  A14.  Displacement  spectrum  of  D event  for  MNV  station 
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Displacement  spectrum  of  E event  for  LAC  station 
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Figure  A17.  Displacement  spectrum  of  E event  for  MNV  station 
( loge  2.7,  Hz  = 0) . 
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Figure  A18.  Displacement  spectrum  of  F event  for  LAC  station 
(loge  2.7,  Hz  = 0) . 


